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ABSTRACT 

We have made a quantitative calculation for the systematic evolution of average extinction by inter- 
stellar dust in host galaxies of high-redshift Type la supernovae, by using a realistic model of photometric 
and chemical evolution of galaxies and supernova rate histories in various galaxy types. Wc find that 
average B band extinction (Ab) at z '^ 0.5 is typically 0.1-0.2 mag larger than present, under a natural 
assumption that dust optical depth is proportional to gas column density and gas metallicity. This 
systematic evolution causes average reddening with E{B — y) ~ 0.025-0.05 mag with the standard ex- 
tinction curve, and this is comparable with the observational uncertainty of the reddening of high-redshift 
supernovae. Therefore, our result does not contradict the observations showing no significant reddening 
in high-z supernovae. However, the difference in apparent magnitude between an open universe and a 
A-dominated flat universe is only ~ 0.2 mag at z '^0.5, and hence this systematic evolution of extinction 
should be taken into account in a reliable measurement of cosmological parameters. Considering this un- 
certainty, we show that it is difficult to discriminate between an open and A-dominatcd flat cosmologies 
from the current data. 



Subject headings: cosmology: observations 
general 

1. INTRODUCTION 



dust, extinction — galaxies: evolution — supernovae: 



Type la supernovae (SNe la) have been known as a 
representative standard candle in the universe, and used 
in measurements of cosmological parameters such as the 
Hubble constant (Hq), density parameter {^m), and the 
cosmological constant {ft\). Recently two independent 
groups have obtained the same result that a A-dominated 
flat universe is strongly favored, by using a few tens of SNe 
la at redshift z -- 0.5 (Riess et al. 1998 [R98]; Perlmutter 
et al. 1999 [P99]). Extinction by dust could affect these 
analyses, and both groups have made a considerable ef- 
fort to assess the systematic uncertainty due to extinction. 
Both groups reported that there is no significant color dif- 
ference between the high-z and local samples. In the sam- 
ple of P99, average reddening {E{B - V)) is 0.033±0.014 
mag for the local sample and 0.035±0.022 for the high- 
z sample (P99). The mean B — V color of the high-z 
R98 sample is -0.13±0.05 or -0.07±0.05 depending on 
two analysis methods, while expectation of unreddened 
color is —0.10 to —0.05. However, there is a statistical 
uncertainty of > 0.03 mag in the color difference, and 
in addition, there is a systematic uncertainty in the K- 
correction of about 0.03 mag (P99, and probably similar 
number also for R98). Therefore, there is observational 
uncertainty of at least ^ 0.05 mag in reddening evolution 
for both groups. Furthermore, these error estimates have 
been achieved by statistical averaging of many supernovae; 
typical color uncertainty for each supernova is much larger 
(~ 0.1-0.2 mag, see Fig. 6 of P99). Therefore, it is difficult 



to clearly check a systematic evolution of average B-band 
extinction AAb ~ RbAE{B - V) ^ 0.2 mag with a red- 
dening of AE{B — V)^ 0.05 by the current observational 
data, where Rb is the total-to-selective extinction ratio for 
the B band, and Rb ^ 4 for the standard extinction curve. 
On the other hand, the difference of apparent magnitude 
between an open universe and a A-dominated flat universe 
is only ^ 0.2 mag at z ~ 0.5, and hence it is still possi- 
ble that unchecked systematic evolution of extinction has 
affected the measurements of cosmological parameters. 

Aguirre (1999a, b) considered the effect of intergalac- 
tic dust which was ejected from galaxies. Such dust could 
have a significantly grey extinction and Aguirre has shown 
that this kind of dust may affect measurements of cos- 
mological parameters with smaller reddening. This is an 
interesting possibility, but existence of such intergalactic 
dust has not yet been confirmed. Here we consider the 
dust existing in host galaxies with the standard extinction 
curve. Although it is uncertain whether supernovae evolve 
to high rcdshifts, their host galaxies should undoubtedly 
evolve as shown by various observations of galaxies at high 
rcdshifts. Both the gas column density and gas metallicity, 
which are important physical quantities for dust opacity, 
change with time by various star formation histories de- 
pending on morphological types of galaxies. Therefore, 
average dust extinction in host galaxies should evolve sys- 
tematically, and the aim of this paper is to make a quan- 
titative estimate for this evolution, by using a realistic 
model of photometric and chemical evolution of galaxies 
and supernova rate histories in various types of galaxies. 



We find that typical evolution in average Ab is ^ 0.1-0.2 
mag from z — to ^^0.5, which is significant for nieasure- 



ments of cosmological parameters but may have escaped 
from the reddening check. Therefore, this effect should not 
be ignored in measurements of cosmological parameters by 
high-z SNe la. 

2. EVOLUTION OF AVERAGE EXTINCTION IN 
HOST GALAXIES 

In this letter we consider only the average extinction 
of a supernova in a host galaxy, and do not consider the 
variation within a galaxy depending on the supernova lo- 
cation in it. Although the variation within a host galaxy 
can be washed out by statistical averaging of many super- 
novae, evolution of galaxies will cause systematic evolution 
of average extinction which cannot be removed by statis- 
tical averaging. It is physically natural to assume that the 



formation stops at about 1 Gyr after the formation by a 
supernova-d riven galactic wind (Larson 1974; Arimoto & 
Yoshii 1987). We assume that gas fraction in a ellipti- 



cal galaxy decreases exponentially after the galactic wind 
time {r^ 1 Gyr), with a time scale same as the galactic 
wind time. These models give the evolution of gas frac- 
tion and metallicity in each galaxy type depending on the 
star formation history. 

SN la rate history in each type of galaxies is calculated 
with the metallicity-dependent SN la model introduced by 
Kobayashi et al. (1998). In their SN la progenitor model, 
an accreting white dwarf (WD) blows a strong wind to 
reach the Chandrasekhar mass limit. If the iron abun- 
dance of progenitors is as low as [Fe/H] ^ — 1, the wind 
is too weak for SNe la to occur. Their SN la scenario 
has two progenitor systems; one is a red-giant (RG) com- 
panion with the initial mass of M^q.a '^ ^Mp, and an 



dust-tj O - mctal ratio is constant and hence the duat opacity 
is pro portional to gac column density and gae metallicity 



of a host galaxy. In fact, it is well known that the extinc- 



orbital period of tens to hundreds days (Hachisu, Kato, & 
Nomoto 1996| , 1999). The other is a near main-sequence 
(MS) companion with an initial mass of Mms,o ~ 2-3Mp 



tion ii n our Galaxy is well correlated to the III gas column 
densitl }^ (e.g., Burotoin fe Hoiloo 1982; Poi 1902). It io aloo 



known that the dust opacity is correlated to the metallic- 
ity among the Galaxy and the Large and Small Magellanic 
Clouds, when gas column density is fixed (e.g., Pei 1992). 
Hence in the following we assume that the dust opacity 
is proportional to gas column density and gas metallicity, 
which evolve according to the star formation history in a 
galaxy. 

The star formation history can be inferred from the 
present-day properties of observed galaxies, by using the 
well-known technique of stellar population synthesis. We 
can estimate the time evolution of gas fraction and metal- 
licity in a galaxy by using photometric and chemical evo- 
lution models for various galaxy types. Since the observed 
extinction of high-z SNe la is an average over various types 
of galaxies, we also need the evolution of SN la rate in 
various galaxy types. In the next section, we describe the 
model of galaxy evolution and SN la rate evolution used 
in this letter, which is constructed to reproduce various 
properties of the present-day galaxies. 



2.1. 



Evolution of galaxies and Type la 
supernova rate 



We use photometric and chemical evolution models for 
five morphological types of E/SO, SOa-Sa, Sab-Sb, Sbc-Sc, 
and Scd-Sd. The basic framework of the model is the same 
as that of elliptical galaxies of Arimoto & Yoshii (1987) 
and that of spiral galaxies of Arimoto, Yoshii, & Taka- 
hara (1992), but model parameters are updated to match 
the latest observations (Kobayashi et al. 1999), by using 
an updated stellar population database of Kodama & Ari- 
moto (1997) and nucleosynthesis yields of supernovae of 
Tsujimoto et al. (1995). The model parameters for spi- 
ral galaxies are determined to reproduce the present-day 
gas fractions and B — V colors in various galaxy types 
at 15 Gyr after the formation. The model of elliptical 
galaxies is the so-called galactic wind model, in which star 

^In reality, there should be some dispersion in galaxy a^ 
ing to the fact that all galaxies should become younger on 
cannot remove this systematic effect. 



and a period of several t enths of a day to sev eral days 
& van den Heuvel 1997; Hachisu et al. 1999). The occur- 
rence of SNe la is determined from two factors: lifetime 
of companions (i.e., mass of companions) and iron abun- 
dance of progenitors. (See Kobayashi et al. 1998, 1999 for 
detail.) This model successfully reproduces the observed 
chemical evolution in the solar neighborhood such as the 
evolution of oxygen to iron rat io and the abundance d istri- 
bution function of disk stars ( pK^obayashi et al. 1998[ ), the 
present SN II and la rates in spirals and ellipticals , and 
cosmic SN la rate at z ^ 0.5 ( [Kobayashi et al. 1999| ). 

2.2. Average extinction evolution towards high 

redshifts 

We have modeled the evolution of gas fraction (/g), 
metallicity (Z), and SN la rate per unit baryon mass of 
a galaxy {TZja) in various types of galaxies, from which 
we calculate the evolution of average extinction in the 
universe. We assume that these quantities do not de- 
pend on the mass of galaxies. The basic assumption is 
that the dust opacity, and hence average Ab in a galaxy 
is proportional to gas column density and gas metallic- 
ity. The average extinction at redshift z in a i-th type 
galaxy with the present-day B luminosity Lb is given 
by AbAz, Lb) = ^iU^{t^)Z^{t^)[r^4LB)]-^{^'h/LB)^LB, 
where t^ is the time from formation of galaxies, re the ef- 
fective radius of galaxies, and {M^/Lb) is the baryon-mass 
to light ratio which is determined by the evolution model. 
We assume a single formation epoch zp for all galaxy types 
for simplicity.^ The proportional constant k will be de- 
termined later. We do not consider the size evolution of 
galaxies, and determine re{LB) from empirical relations 
observed in local galaxies (Bender et al. 1992 for ellip- 
ticals, and Mao & Mo 1998 for disk galaxies). It should 
be noted that the extinction depends on the absolute lu- 
minosity of galaxies. From the empirical Lb-Vc relation, 
the surface brightness becomes brighter with increasing lu- 
minosity of disk galaxies, and hence the massive galaxies 

;es. However, the systematic evolution of dust extinction is ow- 
average towards high redshifts, and present-day age dispersion 



should be more dusty than smaUer ones. This trend is 
consistent with observations (van den Bergh and Pierce 
1990; Wang 1991). Then the average extinction of SNe la 
over all galaxy types at a given redshift is 



{Ab{z)) 



J2J dLBAB4z,LB)nl,4tz)iMh/LB)^LBMLB) 
Y.^ ! dLBnl,,^{tz){M^/LB)^LB(t>^{LB) 



where (pi is the type-dependent galaxy luminosity function 
at z = 0, for which we adopted the Schechter parameters 
derived by Efstathiou, Ellis, & Peterson (1988) using the 
catalog of the Center for Astrophysics (CfA) Redshift Sur- 
vey (Huchra et al. 1983). 

We have to determine the overall normalization of ex- 
tinction, K, for which we use the average V extinction of 
the Milky Way (Sbc type, Lb = 1.4 x IO^^Lbq) at z = 0: 
{Av)mw- This is an average of extinction of SNe la occur- 
ring in our Galaxy seen by an extragalactic observer, and 
hence it is different, in a strict sense, from the average of 
the Galactic extinction which is extinction of extragalac- 
tic objects observed by us. However, if the location of 
the Sun is typical in the Milky Way, we may infer this 
quantity by the average of the Galactic extinction. The 
average Galactic extinction of the 42 SNe la observed by 
P99 is ~ 0.1 mag in Ar or Aj (see Table 1 of P99). This 
suggests {Ay)mw ^ 0.1-0.2 with the standard Galactic 
extinction law (e.g., Pei 1992). The average reddening 
of the Galaxy then becomes {E{B - V^))mw ~ 0.03-0.06 
mag, which is a typical reddening at the Galactic latitude 
of ^40-50° in the Galactic extinction map (Burstein & 
Heiles 1982; Schlegel, Finkbeiner, & Davis 1998). This es- 
timate is consistent with a model of dust distribution in 
our Galaxy, which suggests that the average extinction of 
SNe la in the Galaxy seen by an extragalactic observer 
is typically {Av)mw ~ 0.1-0.2 mag (Hatano, Branch, 
& Deaton 1998). ^ Therefore we use {Av)mw ~ 0.1-0.2 
mag as a plausible range of the average extinction of our 
Galaxy. 

Figure 1 shows the evolution of B extinction for each 
galaxy type as well as the average over all galaxy types, 
normahzed by {Av)mw, i-e., {Ab{z)) / {Av)mw ■ Here we 
used a cosmological model with {h, V,m, ^a) = (0.5, 0.2, 0), 
and set zp = 4.5 so that the age of galaxies is 15 Gyr which 
was assumed in the evolution model. The thick solid line is 
the average over all galaxy types, and the thin lines are for 
individual galaxy types as indicated. Since we have used a 
galactic wind model for elliptical galaxies, they do not have 
interstellar gas and hence there is no extinction in elliptical 
galaxies at z < 1 . The evolution of extinction is caused by 
spiral galaxies, but the behavior of evolution is consider- 
ably dependent on galaxy types. Early-type spiral galaxies 
become more dusty towards z ~ 1, but an opposite trend 
is seen for late types. These behaviors can be understood 
as a competition of the two effects: gas fraction evolution 
and metallicity evolution. The gas fraction in early spiral 
galaxies is much smaller than late types at present, but 
rapidly increases towards high redshifts. This increase is 
responsible for increase of gas column density and hence 

^Hatano et al. suggested that most supernovae are only mildly obscured but there is a long tail to stronger extinctions 
in the extinction distribution. In the actual observations, such a tail will be cut out due to a magnitude limit of a survey. 
Hence, we have used here the mean extinction of the "extinction-limited subset" in the Table 1 of Hatano et al., in which 
strongly obscured supernovae with Ab > 0.6 are removed. 



the dust opacity. On the other hand, the gas fraction does 
not increase so much in late type galaxies, and decrease 
of metallicity towards high redshifts is responsible for the 
decrease of dust opacity. In redshifts more than 1, the 
extinction decreases towards higher redshifts in all spiral 
galaxies because the metallicity evolution becomes domi- 
nant. 

The average over all types is weighted by the SNe la rate 
in each type. Because the star formation rate increases 
more rapidly to high redshifts in early-type spiral galaxies 
than late types, the average extinction is more weighted 
to early types at higher redshifts. Hence {Ab) / {Av)mw 
increases to high redshifts by ~ 1 from z = to 0.5. This 
result suggests that, with {Av)mw ~ 0.1-0.2, the average 
extinction (Ab) of SNe la at z ~ 0.5 is larger than the local 
sample by about 0.1-0.2 mag. This systematic evolution 
of average extinction is comparable with the difference be- 
tween an open and a A-dominated universe in the Hubble 
diagram, and hence this effect significantly affects mea- 
surements of cosmological parameters. In the next section 
we apply the above model in the estimate of cosmological 
parameters by using the sample of P99. 

3. EFFECT ON THE COSMOLOGICAL 
PARAMETERS 

Figure 2 shows the Hubble diagram for SNe la of 
the primary fit C of P99, which plots restframe B 
magnitude residuals from a A-dominated flat cosmology 
[{h,nM,^A) = (0.65,0.2,0.8)] without dust effect (thin 
solid line). Thin long- and short-dashed lines are the 
predictions of the dust-free case with an open universe 
(0.5, 0.2, 0.0) and the Einstein-de Sitter (EdS) universe 
(0.5, 1.0, 0.0), respectively. As reported by P99, the 
A-dominated flat universe gives the best-fit to the data. 
Next, the thick lines show the predictions when the model 
of extinction evolution is taken into account, where the 
cosmological parameters are the same with the dust-free 
curves of the same line-markings. Here we adopt {Ay)mw 
— 0.2 mag. In the open and A-dominated models, the 
galaxy formation epoch is set to zp = 4.5 and 5.0 so that 
the age of galaxies becomes 15 Gyr. In the EdS model, 
the age of the universe is shorter than 15 Gyr, and hence 
we set Zi^ = 5 which gives an age of galaxies as ~ 12 Gyr. 
Although this age is a little shorter than that assumed in 
the evolution model, the evolutionary effect during 12-15 
Gyr is small and hence this inconsistency is not serious. 
As expected, the model curves with the dust effect are 
typically 0.1-0.2 mag fainter than those without dust. As 
a result, the open universe becomes the most favored cos- 
mology among the three when the extinction evolution is 
taken into account. 

We avoid more detailed statistical analysis to derive any 
decisive conclusion about the cosmological parameters, be- 
cause the result would be highly dependent on the extinc- 
tion evolution model. However, the evolution model pre- 
sented here is quite a natural and standard one without 
any exotic assumption. Hence, our conclusion is that the 
systematic evolution of average extinction in host galaxies 



should more carefully be taken into account when one uses 
SNe la to constrain the cosmological parameters. 

4. DISCUSSION 

P99 estimated the systematic uncertainty of extinction 
to be less than 0.025 mag (Ict) in Ab, based on an analysis 
after removing nine reddest supernovae. Aguirre (1999b) 
argued that this limit does not apply if the dispersion in 
brightness and/or colors of high-z supernovae is dominated 
by factors other than extinction. As noted before, the ob- 
servational uncertainty in E(B — V) for each supernova 
is typically ^ 0.1-0.2 mag, which is larger than the sys- 
tematic reddening of E{B — V) ^ 0.025-0.05 considered 
in this letter. Therefore, it is doubtful that P99 analysis 
successfully removed high-z supernovae reddened by the 
systematic extinction evolution. Rather, the supernovae 
removed by P99 might be reddened by strong reddening 
depending on locations in host galaxies, or by some other 
reasons as argued by Aguirre (1999b). Hence, we consider 
that the systematic evolution of dust considered here has 
not yet been checked by the observations. 

In the analysis of R98, the reddening correction is sys- 
tematically included in the process of light-curve-shape fit- 
ting, and one may consider that our result does not apply 
to the analysis of R98. However, as noted in Introduc- 
tion, one must correct a reddening with E{B — V) ^ 0.05 
to discriminate between an open and A cosmologies, and 
this is comparable with the observational error of color 
estimates. In principle, it is difficult to correct the red- 
dening effect when the reddening is as small as the obser- 
vational error in colors, because the reddening correction 
is performed based on the observed colors. The reddening 
correction by R98 will be effective for supernovae strongly 
reddened beyond the color uncertainty, but it is not clear 
to us whether this correction has successfully corrected the 
systematic reddening evolution discussed in this letter. 



The two groups (R98, P99) argued that the observed 
dispersion of apparent magnitudes showing no significant 
evolution to high redshifts gives further support that their 
results are not affected by extinction. Their argument on 
the dispersion test is true if the observed dispersion is dom- 
inated by the dispersion of extinction, and R98 claimed 
that the observed dispersion is smaller than that expected 
from a dust distribution model of Hatano et al. (1998) 
if mean extinction is ~ 0.2 mag. However, there is large 
uncertainty in dust distribution models, and the selection 
effect might have reduced the apparent dispersion (see the 
footnote 2). We should also be careful about the uncer- 
tainty in the observed 'intrinsic' dispersion, which was de- 
rived by subtracting the measurement error from the actu- 
ally observed dispersion (P99). Hence, this test does not 
favor our scenario but cannot strictly reject it. 

We suggest that the best way to constrain the cosmologi- 
cal parameters by high-z SNe la would be an analysis using 
only SNe la in elliptical galaxies, in which the dust evolu- 
tion effect is much smaller than spiral galaxies at z < 1. 
In fact, P99 tried to analyze their supernovae with known 
host-galaxy types, and found no significant change in the 
best-fit cosmology. However, the host-galaxy classification 
is only based on spectra of host galaxies without high- 
resolution images. The uncertainty of a fit with a specified 
host-galaxy type is still large due to the limited number of 
supernovae, and P99 concluded that this test will need to 
await the host-galaxy classification of the full set of high-z 
supernovae and a larger low-z supernova sample. It is ex- 
pected from our calculation that the fitting result of such 
an analysis in the future will be dependent on host-galaxy 
types, giving important information for chemical evolution 
of galaxies. High-z supernovae beyond z -^ 1 are also de- 
sirable to study galaxy evolution as well as cosmological 
parameters. 
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Fig. 1. — Average _B-band extinction of type la supernovae as a function of redshift. Extinction Ab is normalized 
by the average V extinction of our present-day Galaxy, (Ay)Mw (s66 text). The thick solid line is the average over the 
five morphological types of galaxies, which is weighted by SN la rate in them. Five thin lines are extinction evolution in 
individual galaxy types, as indicated. An open universe with (/i, i^M, ^a) = (0.5, 0.2, 0.0) is assumed, and the formation 
epoch of galaxies is set to zp = 4.5. 
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Fig. 2. — The Hubble diagram of type la supernovae. The data are those used in the primary fit C of Perlmutter 
et al. (1999). Restframe B Magnitude residuals are from a A-dominated universe with (ft,, J^Mj^a) = (0.65,0.2,0.8) 
without dust effect (thin solid Hne). Thin long- and short-dashed Hues are predictions of an open universe (0.5, 0.2, 0.0) 
and Einstein-de Sitter universe (0.5, 1.0, 0.0) without dust. The three thick lines are predictions for the case with dust 
evolution, where the cosmological parameters are the same with the dust-free curves of the same line-markings. 



